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We study cosmology on a conical brane in the six-dimensional Einstein-Maxwell-dilaton system, 
where the extra dimensions are compactified by a magnetic flux. We systematically construct 
exact cosmological solutions using the fact that the system is equivalently described by (6 + n)- 
dimensional pure Einstein-Maxwell theory via dimensional reduction. In particular, we find a power- 
law inflationary solution for a general dilatonic coupling. When the dilatonic coupling is given by 
that of Nishino-Sezgin chiral supergravity, this reduces to the known solution which is not inflating. 
The power-law solution is shown to be the late-time attractor. We also investigate cosmological 
(■ — ' tensor perturbations in this model using the (6+n)-dimensional description. We obtain the separable 

f"*^ , equation of motion and find that there always exist a zero mode, while tachyonic modes are absent 

f^**) 1 in the spectrum. The mass spectrum of Kaluza-Klein modes is obtained numerically. 
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I. INTRODUCTION 



Higher dimensional models have a longstanding history, but there has in particular been a growing interest since the 
JL , discovery of D-branes in string theory, leading to the concept of "braneworld" scenarios. The braneworld framework 
qj- allows us to consider large (or even infinite) extra dimensions due to localization of matter fields on the brane 

The brane models with two extra-dimensions have drawn considerable attention among others. This is because six- 
L - - dimensional (6D) models may be able to eliminate the hierarchy problem in particle phenomenology, and at the same 
time bring us the possibility to detect signatures of extra dimensions, as the compactification radius can be of order 
0.1 mm |2J. Moreover, recent realizations of such brane models with rugby-ball (or football) shaped extra-dimensions 
provide a possible self-tuning mechanism to resolve the cosmological constant problem H , though the mechanism has 
, been criticized for several reasons [3] ■ 

In contrast to the 5D (codimension 1) brane models (e.g., @, @j), cosmological aspects of 6D models have not been 
explored much (e.g., @, §)■ Taking into account the gravitational backreaction of branes, one immediately faces the 
■ problem of the localization of matter: it is difficult for a codimension 2 defect to accommodate energy-momentum 
tensor different from pure tension. This fact hampers attempts to construct 6D models for which 3-branes have 
Friedmann-Robertson- Walker geometry. One will notice, however, that not only brane localized matter fields but also 
those in the bulk can support a cosmic expansion of branes. We focus here on this latter possibility. 1 
^ , The purpose of this paper is to present exact cosmological solutions in the presence of a scalar field, flux, and conical 
3-branes in six dimensions. We will be considering the Einstein-Maxwell-dilaton s yst em, and for a certain choice of 
a parameter our Lagrangian reduces to that of Nishino-Sezgin chiral supergravity [17l Ha]. In Ref. [lj|, Mukohyama 
et al. presented a warped braneworld configuration in a more simplified setu p in Einstein-Maxwell theory. The 
perturbation dynamics and stability issue have been extensively studied in [20l |21[ [22l | . In the context of Nishino- 
Sezgin supergravity, similar warped compactification solutions with Minkowski branes were found in [23l . [24], [25[ and 
de Sitter braneworlds were explored in 26]. While none of these solutions have a nontirivial time-dependence of 
metric functions, Tolley et al. [27| found cosmological scaling solutions in 6D chiral supergravity without showing 
explicitly the metric of the internal space (see also Ref. [28| for a numerical construction of cosmological solutions in 
Nishino-Sezgin supergravity). The solutions given in this paper are exact both for the noncompact four dimensions 
and for the compact internal space. After constructing our cosmological solutions, we shall also discuss the behavior 
of tensor perturbations in the present model. 

We will be taking the approach developed in Ref. [29(, which was originally used to study 5D warped brane 
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models with a scalar field [3Cj, l3l|, [32|j [33J, [3J, l35(. In [29[, instead of solving the relevant 5D field equations, the 
authors employed a (5 + n)-dimensional equivalent description without a scalar field. This approach, with the extra 
n-dimensional space playing effectively the role of the scalar field after dimensional reduction, greatly simplifies the 
analysis of cosmological perturbations as well as constructing dynamical background solutions. Here we use a (6 + n)- 
dimensional Einstein-Maxwell system instead of working directly in the 6D model with a scalar field, and show that 
the same mathematical technique as [29( successfully works. 

The structure of the paper is as follows. In the next section we introduce our dimensional reduction approach and 
demonstrate how one can obtain a Minkowski braneworld using the (6 + n)-dimensional description. In Sec. IIIII we 
construct cosmological solutions in the 6D braneworld. The behavior of tensor perturbations is discussed in Sec. IIVI 
Finally, our conclusions are summarized in Sec. [V] 



II. DIMENSIONAL REDUCTION APPROACH TO 6D BRANEWORLD MODELS 



We are interested in the 6D system described by the action 



-^1 [R[ g ] - d a <pd a <p - 2e"^A) - -e^F ab F ab 



(2.1) 



with two or fewer branes. Here F a b denotes the U (1) field strength. For 7=1 the action (|2.ip coincides with 
the bosonic part of Nishino-Sezgin supergravity [l7l [l8| with some fields set to be zero consistently (and with the 
identification of A — > 2gfM^ 6 y where gi is the U(l) gauge coupling). In this paper we consider more general cases 
and assume that the parameter 7 is arbitrary in the range < 7 < 1. 

The above Einstein-Maxwcll-dilaton system with branes has an equivalent (6 + n)-dimensional description in terms 
of pure Einstein-Maxwell theory. To show this, we start with the (6 + n)-dimensional Einstein-Maxwell action 



j6+n 



X\ 



M, 



4+n 



(TZ [Q] - 2A) - 1 -TmnT 



AIN 



and the (4 + n)-dimensional pure tension branes (i = 1,2) 

s { h t +n) = -f d^x^gVa,, 

where <7t>i is the determinant of the induced metric on the brane. 
Let us consider the metric of the form 



dX M dX N 



-n<t>{*)/2 gab ( x)dx a dx b + e mx) 5mndy ™ dy n^ 



11 D 



and the field strength with 

Tab = Fab{x) and T mU = 0. 
With the above ansatz, dimensional reduction yields 



R[g] 



" ( " 4 +4) daW - 2e-"^ 2 A^) - \e n ^ 2 F ah F ah 



(2.2) 



(2.3) 



(2.4) 



(2.5) 



(2.6) 



where := Af* 6 4 f ™^V( n ), F ao := !F a bV}£)i an( l ^(n) i s the volume of the n-dimensional flat space. Note that in (|2.6 
indices are raised and lowered with g ao . Now, defining 



y/n(n + 4) 
V ■= o <P> 7 : = 



(2.7) 



we see that (|2.6I) is identical to the action (|2.ip . We also define the 3-brane tension as := TiVi n ). Then the brane 
action (12.311 reduces to 



(4) _ 



-Id x^j-qi Ti, 



(2. 
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where is the determinant of the induced metric on the 3-brane. The tension of the 3-brane does not couple to the 
dilatonic scalar field after dimensional reduction. Thus, we may use the (6 + rt)-dimensional pure Einstein-Maxwell 
description (|2.6| with (|2.3p instead of working directly in the 6D Einstcin-Maxwell-dilaton system (]2.1|1 with tensional 
branes (|2.8p . Given a metric Gain of the form (|2.4p and a gauge field that satisfies (|2.5p , one automatically obtains 
gab, Fab, and (p that solve the action (|2.ip . Since the action (|2.2[) is more simple than (|2.ip . it is often easier to find 
a solution of (|2.2j) . Note that in the action (|2.6|) n is an arbitrary parameter and may no longer be an integer. The 
right signature for the scalar kinetic term [n(n + 4) > 0] ensures that j 2 is nonnegative. 
The field equations derived directly from the action (|2.ip possess the scaling symmetry 



9 a b 



ip — ► ip + 7 1 lnw, 



with constant u. Note that this is not a symmetry of the action, as 
point of view, this symmetry is manifest because we may write 



Qmn dX M dX N 



— n(jy/2 — \n u 



ug a bdx a dx b 



(2.9) 

2 S'( 6 ). From the (6 + n)-dimensional 



(2.10) 



and use ug a b and <j> + (2/n) lnu instead of g a b and <j). 

We would like to obtain six-dimensional axisymmctric warped compactification solutions with conical branes at 
the points where the geometry pinches off. To this end we make use of the (6 + n)-dimensional generalization of 
the exact solution given in Ref. [l9l |. which is obtained by a double Wick rotation from the higher dimensional 
Reissner-Nordstrom solution. The metric is written as 



1 



g MN dX M dX N = r 2 g af3 dz a dz 13 + — 



2A 



dr 2 



where 



A 



(4+n)D 



C 



fir) 



(2 + n)(3 + n)A 
and g a fj is any metric that solves the (4 



P 2 f(r)d0 2 



Q 2 



(2.11) 
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r 3+n 



2(3 + n)M^ + n n) V {n) AV 



r 6+2n 



n)-dimensional field equations 
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Raf3[g] = 



-^g a f: 



(2.12) 



(2.13) 



The only nonvanishing component of the field strength is given by 

QP i 



2A V /V^Y 



4+n ' 



(2.14) 



The constant f3 can be absorbed in the angular coordinate 9, but we keep it in the metric. 

We assume that f(r) has two positive roots, r + and r_ = ar + with < a < 1. In terms of the new coordinate 
defined by x '■— r / r +, the metric function can be written as / = r+f(x)> where 



fix) 



1 



(4 + n)(5 + n) 



1 - a 



8+2n 



1 



a 3 + n (l - a 5+n ) 1 



1 - a 3+n x 3 +' 
A 



v 3+n 



,6+2n 



{2 + n)(3 + n)Ar 2 



1 - 



/3+n 



v 3+n 



1 - 



X 



3+n 



(2.15) 



The parameter a characterizes warping of the bulk. In the case of a — 1, the bulk is given by a two-dimensional 
sphere and is often referred to as the rugby-ball compactification. 

The rescaling r 2 _g a pdz a dz^ — > g a pdz a dz^ and A/r^_ — > A allows us to write 



G MN dX M dX N = X 2 g a pdz a dz 13 + -L 



d X 2 

fix) 



1 fix)dB 1 



and 



2A X 



4+n '■ 



(2.16) 



(2.17) 
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where 



Q:= 



3 + n 1 - a 5+n 4 + n A 



2 + n A 



5 + nl 



v 3+n 



1/2 



v (3+n)/2 



(2.18) 



With this rescaling, r + in Eq. (|2.15p is absorbed in A and therefore the solution does not have an explicit dependence 
on r + , while g a p still solves Eq. (|2.13p . 

It is clear from Eq. (|2.18[) that A is bounded from above: 



A < A n 



_ (2 + n)(3 + n) 1 - q 5 +» 
' {a) (4 + n)(5 + n) 1 - a 3 +» ' 



(2.19) 



A large value of A reduces the strength of flux, and the flux vanishes if the bound is saturated. Such a bound on the 
brane curvature scale can also be observed in 6D flux compactification models with de Sitter branes [13, [HJ . 

After dimensional reduction n is a parameter of the 6D theory which is directly related to the dilatonic coupling. 
In particular, Nishino-Sezgin supergravity (7 = 1) is reproduced by taking n — > 00. In this limit the metric function 
f(x) is apparently singular. This is due to a bad choice of the coordinate in the extra direction. A regular expression 
is obtained by using £ := ^ 1 + n / 2 anc [ q, := q/ 1 +"/ 2 j in terms of which one can write the metric in (6 + n) dimensions 
as 



4MN 



dX M dX 



N 



= ^ /{2+n) g al3 dz a dz^ + 1 
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2A 
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where 
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•w)(2 + 3w) 
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A2(2 + w) 



1 - 



and $ := ui/3 with u> := (1 + n/2) 1 , Clearly, h is regular for 



1 - 
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■* 00 {to 



£3(l+u>) 



(2.21) 



(6+n) ^ 



Q 



2A ^(8+3n)/(2+n) 



0). The field strength is now given by 

(2.22) 



We assume that has period 2-k. The constant j3 (or (3) controls deficit angles at \ = 1 (£ = 1) an d X = a (£ = <30 
which are given, respectively, by 



<5i = 2tt 



~l/'(l)l 



2tt 



1-f 1^(1)1 



2vr 



f !/'(«)! 



2tt 



1-f 



(2.23) 



Here the prime stands for a derivative with respect to the argument. The conical deficit corresponds to a codimension 
two brane and the tension is given by % /M^\ = Si . Whereas in general there are two branes in the present model, 
one can choose (3 so that one of the deficit angles vanishes. However, for the warped a^l case it is not possible to 
make both of the deficit angles zero. Hence, the warped model includes at least one brane. The deficit angles satisfy 
the constraint 



2tt -S2 
2tt - 5i 



h'(a) 



fc'(l) 



/'(«) 



/'(I) 



(2.24) 



which therefore places the constraint between the brane tensions. 

It is instructive to present here the simplest example: a Minkowski braneworld. The desired solution is generated 
from the (4 + n)-dimensional Minkowski metric: g a pdz a dz 13 = r\ ilv dx ll dx v + 5 mn dy m dy n (and hence A = 0). We 
identify — \ — ' 2+ ™- 1 and then the 6D metric is 



g ab dx a dx b = e j[lW h^dx v + ^ 



de_ 

HO 



(2.25) 
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The field strength and the scalar field is given, respectively, by 



M (<5)P Q 

F?e = ^2k ^(2+ 7 2)/(l+ 7 2) ' ( 2 - 26 ) 



<P = T-p—2 ln ^- ( 2 - 27 ) 

Since = Ti/M^y the conical deficit and the 3-brane tension are related by Tj/Mfgs = Si. One can check that 

the above configuration reduces to the known solution in the Einstein-Maxwell model [20, [21| and in Nishino-Sezgin 
supergravity [23|, [24], [25| for the corresponding value of 7. 



III. EXACT COSMOLOGICAL SOLUTIONS 



Cosmological solutions can be systematically generated from the following Kasner-type metric, 2 

g a pdz a dz fi = -dt 2 + e 2A ^S ij dx i dx j + e 2B{t) 5 mn dx m dx n . 



(3.1) 



3D 



The explicit form of the time-dependent metric functions A(t) and B(t) is derived in detail in Appendix lAl Looking 
at Eq. and the metric ([2~TB]) [or ([2T2H|l ] with ((HHJ), we identify = xe B = £ 2 /( 2 +™)e s . Then the 6D metric is 

found to be 



j ab dx a dx b = f 2 /( 1 +^ 2 ) \-e nB ' 2 dt 2 + e 2A + nB ' 2 ^dx l dxi 
= e /(lW) [-dr 2 + a 2 (T)5 lJ dx*dxi] + ^ 



iB/2 



2A 



dij 2 
W) 



(3 2 h{C)d6 2 



de 

W) 



(3 2 h{£)d9 2 



(3.2) 



where we have introduced the proper time on the brane, r — J e nB / 4 dt, and defined the scale factors of the external 
and internal spaces as a(r) := e A + nB / 4 an d 5( r ) : — e nB ^ 4 , respectively. The dilaton is given by 



^(T,£) = -ln&(T) + _^lne 

7 1 + 7 

The field strength F^g is the same as that in the static solution ()2.26|) . 

Let us first consider the simple case where the metric functions are given by 

A(t) = H t + A 0l B(t) = H t + B , 

with H% := 2A/[(3 + n)(2 + n)] 7^ 0. In this case it is easy to see that 

a(r) oc r 1 / 7 , b(r) oc t, 



(3.3) 



(3.4) 



(3.5) 



leading to power-law inflation for 7 < 1. However, the Nishino-Sezgin supergravity model, 7=1, does not give an 
accelerated expansion. The 7=1 solution was first found by Tolley et al. [2jj. While they have not presented an 
explicit form of the metric of the two-dimensional internal space, Eq. (|3.2[) shows the full spacetime metric in six 
dimensions. 

In addition to the above trivial solution, we have another type of solutions, 



JA+nB = c inh 



3 + n 
2 + n 



2A (f - 1 ) 



-A-B 



C4 < tanh 



3 + n / A , 



-1 ~> ± V / (2+n)/3r, 



(3.6) 
(3.7) 



2 We derive the main results in this section assuming n > and hence 7^0. Obviously, the n = case does not admit the nontrivial 
Kasner-type solution and only the de Sitter geometry (as in 19]) is allowed on the brane. 



6 



At early times, t — to ~ 0, we have a ~ t p± and b ~ r 9± where 



p± = WT^f) ' q± = 6T5? • (3 - 8) 

The ranges of p± and q ± are 1/3 < p+ < (9 + 4a/3)/33 « 0.48, 1/3 > p_ > (9 - 4\/3)/33 « 0.063, > g+ > 
(1 - 2\/3)/ll « -0.22, and < g_ < (1 + 2\/3)/ll « 0.41. Note that g + is negative and hence in this case the 
compact two-dimensional space is initially contracting. At late times, t — * 00, we have a ~ 7- 1 / 7 and & ~ r. Therefore, 
the solution (|3.5[) is the late-time attractor. This attractor behavior on the 4D brane emerges as a consequence of 
Wald's cosmic no hair theorem [38j for the (4 + rt)-dimensional metric g a p'- the initially expanding Bianchi-type 
metric (|3.ip isotropizes in the presence of the positive cosmological constant A. 

The cosmological dynamics can also be understood in terms of the 4D effective theory in the zero-mode sector. 
In Appendix [B] we show that the induced metric on the brane can be derived from the Einstein-scalar field action 
with an exponential potential up to conformal transformations. Given generic initial conditions for the "position" 
and "velocity" of the scalar field, the scale factor and the time-dependence of ip are obtained from (|3.6p and (|3.7|) . In 
this case the two of the three integration constants correspond to the initial position and velocity of the scalar field, 
whereas the third integration constant corresponds to the overall normalization of the scale factor. The solution (|3.5p is 
realized for the special choice of the initial velocity (depending on the initial position) , as f|3.4|) has only two integration 
constants. The power-law solution obtained from (|3.4p is the late-time attractor for the exponential potential of the 
scalar field. The cosmic expansion looks very different in a different frame, but the discussion here in terms of the 4D 
effective theory (|B3[) is still helpful. 

The dynamics of the A — model is the same as the early-time behavior of the general A ^ case: a ~ t p± and 
5 T i± . From Eq. (|3.8p we see that the 7 = 1 case gives the scale factor identical to that in [36|, [37| ■ One should 
note here that the geometry on the brane is independent of the warping parameter a and the brane tension. They 
only affect the bulk geometry. For this reason, it is not surprising that we can reproduce the scale factor of the model 
without warping and branes [36l |37| just by taking 7 = 1. 



IV. COSMOLOGICAL TENSOR PERTURBATIONS 

We will be discussing the behavior of perturbations in the cosmological background solution constructed in the 
previous section. Here we make use of the (6 + rt)-dimensional description again, and focus on the axisymmetric 
tensor perturbations. 3 The perturbed metric is given by 



mn + Wmn) dX M dX N = x 2 [~dt 2 + e 2A + fty) dxW + e 2B 6 mn dy m dy n ] + -L 



dx 



2 



(4.1) 



2A 

where hij = hij (t, x, x) and = d J hij = 0. The (6 + n)-dimensional Einstein equations give 

e ~3A-nBQ t ^ 3A +nB dth .^ _ ^A^.. = 2 A X - 2 - n 8 x ( X 4+ "/5 x %) , (4.2) 

which is separable. In the Fourier space, the general solution can be written as 

hij = ]>> m WO„(x)e£ T (k;x), (4.3) 

m 

where eJ, T is the transverse traceless tensor harmonics which satisfies V 2 e^ T = — k 2 ej^- . 
The mode solutions ip m and Vt m obey 



d 2 („ : ^ d k 2 2 



d* + [ 3A + nB )dt + ^ +fl _ 
1 d ( A+n t d \ ( fi 2 



i> m = 0, (4.4) 



x ^f^-n m + = o, (4.5) 



X 2+n d X V d X m J 2A 



Perturbations are decomposed into scalar, vector, and tensor parts according to their transformation properties with respect to the 
three-dimensional space. Scalar and vector perturbations will be investigated in a forthcoming paper. 
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where the dot denotes a derivative with respect to t and fi m is a separation eigenvalue. At the poles the metric 
function / vanishes. Therefore, the boundary conditions for Q m are 



2 t 'rV i f^m o 







at 



X = 1, a- 

Two things should be remarked here. First, there always exist a zero mode (/!q = 0): 

Oo = const. 



(4.6) 



(4.7) 



Second, there are no tachyonic modes (fi m < 0) in the spectrum. This can be shown as follows. Using Eq. (I4.5|) we 
obtain 



2A 



x 2+n n 2 m d x 



o ro ( x 4+n /n^)' = / x i+n f (Kf d x > o 



(4.8) 



where we used /(l) = f{a) = and assumed the regularity of f2 m and fl' m at the boundaries. Thus, fi m is not 
negative. The equality holds only for Q' m = 0, which is definitely the zero mode. 
Eq. (|4.4[) can be written in terms of the conformal time coordinate r) := J e~ A dt as 



dr\ 2 



ndB\ d 
2 dr] J dij 



^+[2H+^)^ + k 2 + e- nB / 2 ^ m a 2 



4>m = 0, 



(4.9) 



where H := a 1 da/dr] is the conformal Hubble parameter on the brane. We see from Eq. (|4.9[) that the Kaluza-Klein 
(KK) mass with respect to observers on the 4D brane is time-dependent and is given by 



-7iB(7))/2 2 



(4-10) 



As for the zero mode, Eq. (|4.9[) gives the standard 4D equation for a massless graviton except for the extra term 
dB/dij. The presence of this term reflects the fact that the effective theory for the zero-mode sector is given by a 
scalar-tensor gravity, as is argued in Appendix [Bl 

We are mainly interested in the inflationary attractor solution (|3.4[) . In this case, Eq. (|4.9[) reduces to 



2 + n d 



Ipm = 0, 



^drf 2 T] dr) H^rj 2 t 

which can be solved exactly for any value of /i^ n . The general solution is given by 

i> m = (-M (3+n)/2 ^^(-M, 

where Z,, is a linear combination of Bessel functions of order v and 



H 2 



(3 + nf 



(4.11) 



(4.12) 



(4.13) 



The late-time asymptotic behavior of (|4.12[) depends on the value of p? m . "Heavy" modes with fj, m > fi 2 , where 

(3 + n)A 



2(2 + n)' 



(4.14) 



decay rapidly as \ip m \ ~ a 2 ( 3 +")/( 4 +") at late times (— kr\ — > 0). "Light" modes with ^ < /i 2 decay more slowly, 
and for p? m <C ^ we nave IV'ml ~ a^m/^e '( 3 +n)/( 4 +™). j n particular, the zero-mode fluctuations are frozen, ipo ~ 
const., at late times. Since 



Mia 2 

n 2 



2(3 + n) 



"1 2 ,,2 



(4 + n) 



2 ' 



(4.15) 



heavy (light) modes in the (6 + n)-dimensional description are also heavy (light) for a 4D brane observer. 



a 



a 



FIG. 1: Eigenvalues of the first three Kaluza-Klein modes, /if, f/%, and /j§, versus a for n = 10. The parameter is given by 
A/A = 0.156 (left) and A/A = 0.74256 (right). The dashed line indicates fi 2 c . 

Now we are to determine the mass spectrum of KK modes {/i 2 ,/i2, • ' '} numerically. The spectrum depends on n 
and the parameters of the background solution, a and A. Given a set of parameters, we numerically solve Eq. (|4.5p 
supplemented with the condition ft(l) = 1 and fi'(l) = -/i 2 /[2A/'(l)]. Then we have 



as a function of /z 2 . The mass spectrum is given by zeros of I{^ 2 )- 

Examples of /x^ are shown in Fig. [TJ As is clear, /x^ is an increasing function of a. One finds that /xf can be 
smaller than /x 2 for a relatively large value of A. To see this more closely, we plot in Fig. [2] the lightest KK mass /xf in 
the limit of a — > as a function of A. We find that for n~lwe have /xf ~ /x 2 when A is close to its maximum value, 
while larger n tends to give a very light mode for large A. 

If we have very light modes in the spectrum, such modes can be excited during the period of inflation and decay 
slowly. This is certainly an interesting situation, but is not favored for the following reason. According to Ref. [39j |. a 
large expansion rate in the external coordinates generally has the effect of destablizing de Sitter compactifications. In 
fact, the presence of A (corresponding to the (4 + n)-dimensional Hubble rate) decreases the strength of flux as seen in 



Eq. (|2.18p . It was also confirmed directly in [22J that in the 6D flux compactification model there appears a tachyonic 
mode in the scalar sector of perturbations when the expansion on the brane is too large (but still less than the upper 
bound given by Eq. (|2. 10|1 with n = 0). 4 From this we may expect that the (6 + n)-dimensional Einstein-Maxwell 
model will also be unstable for large A, leading to the instability of the 6D Einstein-Maxwell-dilaton model since the 
two descriptions are equivalent. In light of this, it is natural to consider that smaller A is favored also in the present 
model. For sufficiently small A all the KK modes are heavy enough to decay rapidly at late times. 

For understanding the behavior of KK perturbations in more detail, it would be helpful to provide an analytic 
argument on the mode functions and eigenvalues in the special case of a = 1. Rewriting Eq. (|4.5[) in terms of the new 
coordinate w := (2\ — 1 — aO/(l — °) an( l then taking the limit a — > 1, we get the equation 



/(/x 2 ) : 



a 2 f'Sl' + ^-n 

J OA 



(4.16) 





(4.17) 



subject to the boundary conditions 




at w 



±1, 



(4.18) 



where A := A/A max (l). The mode solution is given up to the overall normalization by 




(4.19) 



4 We should remark that in [2211 scalar, vector, and tensor perturbations are defined with respect to the 4D de Sitter spacetime 
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FIG. 2: First Kaluza-Klien eigenvalue in the limit of a — > versus A/A m ax(0). The dashed line indicates no- 
where P m is the Legendre function of the first kind of order m and m = 0, 1, 2, • • • . We checked that this analytic 
result agrees with our numerical calculations in the a — ► 1 limit. 
It follows from the mass spectrum (|4.19|) that 

4 = 4m(m+l)(l±^i-lY (4.20) 

Very light modes emerge if and only if A w A max (l) and n is very large. 

Finally, we shall make a brief comment on scalar perturbations in the present model. Our (6 + n)-dimensional 
description would be especially powerful when studying scalar perturbations, as is the case in the Randall-Sundrum- 
type braneworld inflation model with a bulk scalar field 29]. This is mainly because the extra scalar degree of 
freedom in six dimensions is encoded in a metric perturbation in (6 + n)-dimensions. This fact allows one to find a 
convenient gauge in which both perturbation equations and boundary conditions become simple enough to handle. 
The emergence of a tachyonic perturbation mode in the scalar sector in the 6D Einstein-Maxwell model [12] implies 
that such instability may arise in general (6 + n)-dimensions, and hence the scalar perturbations need a detailed 
investigation (see also Refs. [13, Sl[ for scalar perturbations in the Minkowski brane models in the context of Nishino- 
Sezgin supergravity). This issue will be considered carefully in a forthcoming publication. 



V. SUMMARY 



In this paper we have proposed a new systematic method quite useful for studying brane cosmology in warped flux 
compactifications. The crucial key is that the Einstein-Maxwell-dilaton system in six dimensions is equivalent to the 
(6 + n)-dimensional pure Einstein-Maxwell system under our metric ansatz. Using this fact, we have demonstrated 
that time-dependent warped compactification solutions can be obtained easily. Our "seed" is essentially the double 
Wick rotated black hole solution in Einstein-Maxwell theory, which we believe is more simple and familiar. For 
a general dilatonic coupling 7, we have found a power-law inflationary solution and two nontrivial time-dependent 
solutions. The former solution was turned out to be the late-time attractor for A ^ 0. In the case of Nishino-Sezgin 
supergravity (7 = 1), the power-law solution reduces to the noninflating solution of [13] • Here we have obtained the 
explicit form of the metric of the internal two-dimensional space. For the particular choice of the parameter, A = 0, 
we have seen that our solution reproduces the cosmological solution of [36l . [37j (but with warped bulk geometry and 
conical deficits, which do not affect the geometry of the external space). 

We have also investigated the dynamics of tensor perturbations using our (6 + n)-dimensional description. We 
obtained the separable equation of motion for axisymmetric tensor perturbations, and showed that there always 
exists a zero mode, whereas we have no tachyonic modes in the spectrum. For relatively small A, which is likely from 
the stability consideration, KK modes are too heavy to be excited during inflation. 
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APPENDIX A: DERIVATION OF THE KASNER-TYPE METRIC 



In this appendix we present a detailed derivation of the Kanser-type metric used in the main text. The metric takes 
the form 



g afl dz a dz 13 = -dt 2 + e 2A{t) 5 lj dx l dx ] +e 2B ^8 mn dx m dx n . 



(Al) 



The field equations R a p = [2/(2 + ri)]Ag a p read 

(tt): 3 + i 2 ) +n(B + B 2 ) 

(ij): A + A (sA + nB 

(ran): B + B (zA + nB 



2A 
2 + n 
2A 

2 + n 

2A 
2 + n 



(A2) 
(A3) 
(A4) 



where the dot denotes a derivative with respect to t. From Eqs. (|A3[) and (|A4[) we get ( e 3A + nB \ = \e 3A+nB and 



(A- 6) 



,3A+nB 



= 0, where A := [2(3 + n)/(2 + n)]X. For A 7^ 0, integration of these two equations gives 

e 3A+nB = Cie Vlt + c 2 e^\ (A5) 

-B)e^ B =C 3] j^±^. (A6) 
Substituting this result into Eq. (|A2|) . we obtain 



4CiC 2 + Cl = 0. 



(A7) 



If C3 = 0, the relation (|A7|) implies that C2 = (We are interested in an expanding case). Then we have the de Sitter 
solution: 

A(t) =H„t + A , B{t)=H t + B , (A8) 

with Hq = 2A/[(3 + n)(2 + n)]. This solution has two integration constants, Ao and Bo. If C3 7^ 0, we may write 
Ci = ±(C 3 /2)e~^*° and C 2 = T(C 3 /2)e^ to , where t is a constant. In this case we have 



„3A+nB _ 



C3 sinh 



z A - B = Ca I tanh 



X(t-t ) 



(t - t Q ) 



±^(2+n)/3n 



(A9) 
(A10) 



where we have introduced another integration constant C4 in addition to C3 and to- This solution has three integration 
constants, in contrast to the de Sitter case. 
In the case of A = it is easy to find 



ZA+nB 



= Ci(t-to), 



e A - B = C 2 (t-tof^ {2+n)/3n . 
The behavior of the solution is the same as that of (|A9[) and (|A10I) at t ~ to. 



(All) 
(A12) 



APPENDIX B: EFFECTIVE THEORY IN THE ZERO MODE SECTOR 

The effective four-dimensional action governing the induced geometry on the brane at £ = 1 is 



S e s = 



M (4) 



d 4 Xy/—q 



e nB/2 R (4) _ UBD lL e »B/3qi»> dtt BB v B - 2A 



(Bl) 
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with wbd := ^r 2 -- This can be obtained from the Einstein-Maxwell action (|2.2[) simply by substituting the ansatz 



g MN dx M x N = x 



-nB( X )/2 q ^ x)d ^ dx u + e mx) §mn dy™dy n 



4D 



nD 



1 fd X 2 



2A V / 



j3 2 fd9 2 



(B2) 



and T x q = M?^^j3 /y/2A ■ Q/x 4+n [ E q- PTTTj) ] . with the definition of the four-dimensional Planck mass Mf 



(4) 



The conformal transformation q^ — e nB / 2 q lll , leads to 



M (4) 

Seff = — „ — / d x^-q 



R 



(4) 



«( 2 + ») -uv 



q^ u d^Bd v B - 2Ae 



(B3) 



Thus the zero-mode sector of our model can be described by the Einstein-scalar field system with an exponential 
potential. 
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